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I I In a recent paper [Phys. Rev. E 57, 1550 (1998)] we demonstrated that the 

Sh ' symmetries of the evolution equation and the target state have a profound 

^ ■ effect on the selection of the admissible control parameters. In the present 

I/-) . paper we extend these results to the case of time-periodic target trajectories 

| and inexact symmetries. We also argue that the problem of phase space re- 

\ construction is affected by the presence of symmetries similarly to the control 

!>• ' problem. 
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The desire to improve performance of many practically important systems and devices 
often calls for shifting their operating range into a highly nonlinear area, which after a series 
of bifurcations usually leads to irregular chaotic behavior. This kind of behavior, however, is 
rarely desired, while substantial benefits could be obtained by making the dynamics regular. 
This goal can typically be achieved by applying small preprogrammed perturbations to steer 
the system towards a periodic orbit with desired properties, which is broadly referred to as 
; chaos control. 

The goal of the present paper is to highlight some of the problems that arise in controlling 
^ . systems which possess some kind of symmetry. Even though the importance of symmetries 
in chaotic dynamics has been recognized by a number of authors , symmetric systems 
did not receive adequate treatment in the general framework of chaos control primarily 
because the question of symmetry is largely ignored by the theory of deterministic chaos as 
well as data analysis and control theory. All three disciplines regard symmetric systems as 
nongeneric and, therefore, not very interesting and important. 

However, many practically important dynamical systems, such as spatially extended 
chaotic ones, are symmetric, and thus cannot be successfully treated using the formalism 
developed for generic systems. Indeed, such phenomena as fluid flows, convection or chemical 
reactions often take place inside symmetric containers — cylinders, spheres, pipes and annuli. 
As a result, the dynamical equations also show rotational and translational symmetries. 
Even the dynamics of unbounded systems is often significantly influenced by the symmetries 
of the physical space. Although the presence of symmetries usually simplifies the analysis 
of the dynamics, it also makes system identification and control more complicated due to 
the inherent degeneracies of the evolution operators. In fact, the presence of symmetries, 
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explicit or implicit, makes a number of single-control-parameter methods fail [§]-[§], calling 
for multi-parameter control @J^J^,|]. 

In order to see how the control problem is affected by symmetries, we consider (following 
the analysis conducted in flQ[| ) a general discrete-time system (the arguments for continuous- 
time systems are very similar), whose evolution is described by the map F : I ni xl n " — > M. nx , 
such that 

x' +1 =F(x',u), (1) 

where x* is the n x -dimensional state of the system and u is the n u -dimensional vector of 
system parameters. The objective of control is to make the system follow the (possibly 
unstable) periodic target trajectory x*. Let us linearize equation (H) about this target 
trajectory to obtain 

Ax m = A 1 Ax* + (2) 

where 

A* = D x F(x',u) (3) 

is the Jacobian matrix, which determines the stability properties of the target trajectory, 
and 

B l = D u F(x',u) (4) 

is the control matrix, which defines the linear response of the system to perturbation of 
system parameters. According to linear system theory fLT |, if the target trajectory x* is 



unstable, it can be stabilized by an appropriate feedback through the time-dependent control 
perturbation Au*, provided the matrices A 1 and B t satisfy certain conditions. In the present 
study we concentrate on selecting from the complete set of available system parameters a 
minimal set of control parameters, whose perturbation allows the stabilization of the target 
state, i.e., on making an appropriate choice of the control matrix B l , given the Jacobian 
A 1 . We will see below that the constraints affecting the choice of control parameters can 
be easily obtained from the symmetry properties of the system and the controlled state. 
What is more interesting, symmetry allows one to determine the minimal number of control 
parameters even when the Jacobian A 1 describing the local dynamics is unknown. 

Discrete-time evolution equations of type (P are often obtained as a result of phase space 
reconstruction of a continuous-time system when the dynamical equations 

B(f) = *(8(f),u) (5) 

describing its evolution are unknown. Here s(t) G Q denotes the n s -dimensional inter- 
nal state of the system, and $ is an unknown vector field on the phase space manifold 
Q. Generically, such reconstruction is possible when the measurement of a single scalar 
time-dependent signal y{t), which is a function of the system state s(t), is available. Many 
practically interesting systems, symmetric ones in particular, are, however, extremely non- 
generic and require a number of independent scalar signals for the complete reconstruction. 



Eckmann and Ruelle 12] acknowledged that the choice of signals has to be made carefully 
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by trial and error. Certain general rules concerning this choice, however, can be established 
on purely theoretical grounds, since this problem too can be effectively treated based on the 
knowledge of underlying symmetries []T3| . 



The outline of the paper is as follows. In section [TJ we discuss the implications of sym- 
metry for control of discrete-time linear systems in the vicinity of steady target states. The 
results are generalized for time-periodic target trajectories in section |T|, and for continuous- 
time systems in section |TV]. Section [V] is devoted to the problem of weak symmetry violation. 
The issue of phase space reconstruction is considered in section [VT|. The theoretic results are 
illustrated using several simple examples in section |V11| . Finally, we present our conclusions 
in section |V111| . 



II. TIME-INVARIANT STATES 

A. Stabilizability and Controllability 

Although our analysis is applicable to time-varying systems, we start for simplicity by 
assuming that the target state is time-invariant, x* = x. Then the matrices A 1 and B l 
become constant, and we can drop the time index in @) to obtain 

Ax t+1 = AAx 1 + BAuK (6) 

It is useful to introduce and compare two characterizations of the linearized evolution equa- 
tion (|6]), which extremely simplify the analysis of feedback control algorithms: stabilizability 
and controllability. 

The dynamical system (||) or the pair (A, B) is said to be stabilizable, if there exists a 
state feedback 

Au* = -if Ax*, (7) 

making the system (|6|) stable, i.e., it is possible to find a feedback gain matrix K, such 
that all eigenvalues \' k of the matrix A' = A — BK lie within a unit circle of the complex 
plane, \\' k \ < 1, Wk. Otherwise the system or the pair (A, B) is called unstabilizable. Indeed, 
substituting the feedback (0) into one obtains the linearized evolution equation for the 
closed-loop system 

Ax m = {A- BK)Ax.\ (8) 

with Ax = becoming the stable fixed point of the map (H), if and only if A — BK is stable. 

Since the magnitude of the control perturbation Au* is proportional to the deviation 
Ax* of the system from the target state, feedback of the form (|7D is often called proportional 
in the physics literature, although there are a number of other terms used to denote this 
type of feedback. Control theory uses the term state feedback to refer to the fact that the 
state of the system is used to determine the control perturbation. At first sight equation (|7|) 
seems to impose strict limitations on the allowed form of the feedback law. However, this is 
precisely the form demanded by a number of widely used control algorithms [Q,|],[nj. 

Stabilizability is a property, which usually sensitively depends on the values of system 
parameters. In the majority of practical applications, however, it is preferable to have 
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an adaptive control that would stabilize a given steady state x(u) for arbitrary values of 
system parameters. This is especially important, if one is to track the trajectory x as 
parameters slowly vary, which might be advantageous in many applications, e.g., for moving 
the operating point of a nonlinear device across a bifurcation, from the stable region to the 
chaotic region. Such a control scheme can be obtained, if the more restrictive condition of 
controllability, which is essentially parameter-independent, is imposed on the matrices A 
and B. On the other hand, it can be demonstrated JTT] that the controllability condition 



guarantees that the eigenvalues of the matrix A — BK can be freely assigned (with complex 
ones in conjugate pairs) by an appropriate choice of the matrix K. Therefore, if the system 
is controllable, it is stabilizable as well, and by requiring controllability we satisfy both 
conditions at once. 

The n^-dimensional linear system (0) or the pair (A, B) is said to be controllable if, for 
any initial state Ax' 1 = Axj, times tf—ti > n x , and final state Ax,, there exists a sequence 
of control perturbations Au*% • • • , Au^'" 1 such that the solution of equation (^) satisfies 
Ax'^' = Ax,. Otherwise, the system or the pair (A,B) is called uncontrollable. 

The controllability condition can be represented in a number of different equivalent 
forms. To obtain one particularly convenient form, we make the trivial observation that, if 
it is possible to drive the linear system from an arbitrary initial state Ax,; to an arbitrary 
final state Ax, in n x steps, it is possible to do the same in any number of steps n exceeding 
n x . Suppose we let the system evolve under control for n x steps from the initial state Ax*. 
The final state will be given byQ 

Ax t+n x = (^)n, Ax i + Y^{A) n *- k B&VL t+k -\ (9) 

k=l 

Denote b m the mth column of the matrix B: 

B = [b t b 2 ■•• b n J. (10) 

Regarding the terms (A) nx ~ k h m as vectors in the n x -dimensional tangent space T, 

h k m = (A) n *~ k b m , k = 1, • • • , n x , m = 1, • • • , n u , (11) 

and the control perturbations Ait^* -1 as coordinates, we immediately conclude that equa- 
tion (^]) rewritten as 

Ax, - (A)- Ax, = E (12) 

k=l m=l 

can only be satisfied, if and only if there are n x linearly independent vectors in the set (|TTD, 
i.e., the set {h^} spans the tangent space T. This is equivalent to requiring that 

rank(C) = n x , (13) 



1 Here and below in the text we use the notation (A) n to indicate that A is taken to the power of 
n to differentiate it from the notation A*, where index t defines the time dependence. 
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where the matrix 



C = [B AB (A) 2 B ■■■ {AY^B} (14) 

is called the controllability matrix. Condition ( fL3|) was introduced into the physics literature 
from linear systems theory by Romeiras et al. [| as a simple, but practical test of the 
controllability. 

In contrast, the stabilizability condition requires that the set ( |Tl~D spans only the unstable 
subspace L u C T of the Jacobian A, instead of the whole tangent space T. Stabilizability 
can be formally expressed in the form identical to (f03). Let us define the number of stable 
and unstable^ eigenvalues of the Jacobian n s x and n", respectively (one obviously has n s x + 
n x = n x ). For instance, if A is a diagonalizable matrix, it has n s x linearly independent 
stable eigenvectors which we denote ef, % = 1, • • • , n s x . It can, therefore, be shown using an 
appropriate coordinate transformation that the pair (A, B) is stabilizable if and only if 

rank(iS) = n x , (15) 

where the matrix 

S = [ei ■•• < s B AB •■■ (A) n *~ l B] (16) 

can be called the stabilizability matrix by analogy with the controllability matrix. 

In order to better understand the restrictions imposed on the control scheme by sym- 
metries, it is beneficial to look at the controllability condition from the geometrical point 
of view, assuming n u — 1 and, consequently, B = b. The controllability in this context 
is equivalent to the vectors h 1 , h 2 , • • ■ , h Ux spanning the tangent space T. Generically, the 
matrix A is nondegenerate (has a nondegenerate spectrum) , so one can always find a vector 



b, such that the resulting set (|TT|) forms a basis. However, if A is degenerate, which is a 
usual consequence of symmetry, there will exist an eigenspace of the Jacobian, U C T, 
such that x^A = A r xt , Vx G U with the dimension d r = dim(L r ) > 1, where f denotes 
(complex conjugate) transpose of a matrix or vector. The dynamics of the system in such 
an eigenspace cannot be controlled with just one control parameter (see || for an example 
of such a situation), because the vectors h k only span a one- dimensional subspace of U . 
Indeed, since d r > 1 there will exist d r — 1 adjoint eigenvectors fj G U orthogonal to b and 
each other. Then 

(f,- • h fc ) = i}(A) n *~ k h = A^" fc fjb = A^- fe (f,- • b) = 0, (17) 

so every basis vector h k is orthogonal to every eigenvector f^, j — 1, ■ ■ ■ , d r — 1. 

It is often convenient to define the notion of controllability for individual eigenvectors. 
We will say that the adjoint eigenvector f of the Jacobian A is controllable, if there exists 
m, 1 < m < n u , such that (f • b m ) ^ 0. Respectively, the eigenvector that is orthogonal to 
every column of the control matrix B is called uncontrollable. Using these definitions we 



2 For the purpose of control we regard the central directions, defined by the eigenvalues A such 
that |A| = 1 (Re(A) = in the continuous-time case), as unstable. 
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can, therefore, conclude that the controllability of the linearized system is equivalent to the 
controllability of each and every adjoint eigenvector of the Jacobian matrix (also see [fL5|]). 
Similarly, the stabilizability is equivalent to the controllability of each and every unstable 
adjoint eigenvector. 

If the system dynamics in U happens to be stable (e.g., when the system is stabilizable, 
but uncontrollable), the system can still be stabilized similarly to the nondegenerate case, 
but we have to ensure the controllability in case the dynamics in this eigenspace is unstable. 
This can be achieved by increasing the number of control parameters n u , which extends 
the set (|H|), until it spans every eigenspace of T. This would lead one to assume that the 
minimal value of n u should be defined by the highest degeneracy of the Jacobian matrix A. 
We will see, however, that various kinds of degeneracy have a somewhat different effect on 
the controllability of the system. 

B. Symmetries of the System 

Symmetries usually significantly simplify the analysis of system dynamics, and the control 
problem is no exception. In particular, even when the exact form of the Jacobian matrix 
is unknown, the structure of the symmetry group describing the symmetries of the system 
allows one to reduce the controllability condition ( O ) to a set of much simpler conditions, 
which provide a number of system-independent results. The discussion below is based on 
bifurcation theory and closely parallels the treatment of degeneracy in quantum mechanics 
and spontaneous symmetry breaking in quantum field theory and phase transitions. 

In general we call the system symmetric, if the nonlinear evolution equation preserves 
its form under a set of linear transformations g : x — > x' = g(x) of the phase space. More 
formally, we say that the evolution equation ([!]) possesses a structural symmetry described 
by a symmetry group Q, if the map F commutes with all group actions: 

F(s(x), u) = <7(F(x, u)), Vgeg,VxeT (18) 

or, in other words, if the function F(x, u) is ^-equivariant with respect to its first argument. 
The group Q is usually a byproduct of symmetries of the underlying physical space, such 
as rotational and translational symmetry (domain symmetry), and symmetries of the phase 
space, such as phase symmetry <p <p + 2tt (range symmetry). Since all interesting physical 
symmetries are unitary (such rare exceptions as the Lorentz group are hardly relevant in 
the context of control problem), we will assume that Q is a unitary group. 

Usually, the symmetry demonstrates itself in more than just one way: often steady (as 
well as time-periodic) states x of symmetric systems too will be symmetric with respect to 
transformations g G 7ix, where C Q is an isotropy subgroup of x. In general, the target 
state x might also be symmetric with respect to transformations which do not belong to Q. 
However, considering those does not provide any additional information, so we assume that 

2(x)=x, V<7GW*. (19) 

For the purpose of control it is important to observe that upon linearization about the target 
state x the structural symmetry of the evolution equation ([[]) does not disappear, but is 
replaced with a related dynamical symmetry. Indeed, using the definitions (0), (|T9| ) and 
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the fact that symmetry transformations are linear, one obtains in the linear approximation 
for an arbitrary g G Ti^: 



x + g{AAx) = g(x) + g{AAx) = g(x + AAx.) 

= g(F(x, u) + AAx.) = g(F(x + Ax, u)) 
= F(g(x + Ax), u) = F(g(x) + g(Ax), u) 
= F(x + #(Ax), u) = F(x, u) + Ag(Ax) 

= x + Ag(Ax). (20) 

Defining £ the full symmetry group of the linearized equation @ in the absence of 
control (Au* = 0): 

g(AAx) = Ag(Ax), Wg G C, (21) 

one concludes that the group L describing the dynamical symmetry of the system in the 
vicinity of the target state x includes all transformations g G 7ix, and therefore: 

C C. (22) 

One can speculate that typically C will coincide with H*. As a consequence, if the target 
state x has low symmetry, the symmetry of the evolution equation will be reduced upon 
linearization to a subgroup of Q. However, as we will see in section |V11] , C might be equal 
to Q, or even include Q as a subgroup for highly symmetric target states, with the apparent 
symmetry increased by linearization. 

It turns out that with the help of group representation theory one can substantially sim- 
plify the controllability condition ( |T5j ) and, as a result, obtain a number of useful restrictions 
on the set of control parameters. Consider the matrix representation T generated in the 
tangent space T by the action of transformations g from an arbitrary subgroup CJ of the 
full dynamical symmetry group L: 

(g(yi)) i ^(T(g)x) i ^J2 T ^ j , Vx G T, (23) 

3=1 

where, according to ([H]), all matrices T(g) commute with the Jacobian 

T{g)A = AT{g), Wg G £ C C. (24) 

The knowledge of the representation T is enough to derive a very simple criterion for the 
admissibility of the control matrix. Observe that, if T(g)B = B for an arbitrary transfor- 
mation g G £, then 

C=[T(g)g AT(g)B ■■■ [AY*- l T{g)B] 
= [T(g)B T(g)AB ■■■ T{g){A) n ^~ 1 B] = T(g)C. (25) 

As a result, since Tij(g) ^ 6ij for any g ^ e (where we defined e as the identity transforma- 
tion: e(x) = x), the rows Cj of the controllability matrix become linearly dependent, 

J2( T M - ki)Zj = 0, (26) 
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and the controllability condition fll3|) is violated. Therefore, we obtain a necessary condition 
on the control matrix: 

T(g)B^B, \/geC'\{e}. (27) 

In other words, the control arrangement should be chosen such that the symmetry of the 
linearized evolution equation (|^) is completely broken for (almost all) nonzero control per- 
turbations Au 7^ 0. 



C. Group Coordinates 



Though simple and general, criterion (|27|) is not very helpful for finding the minimal 
set of control parameters satisfying the controllability condition. In order to derive a more 
practically useful criterion one has to make a few more steps. We begin with the reduction 
of the controllability condition to a set of simpler conditions which can be performed [II 



by constructing the Jordan block decomposition of the Jacobian matrix. This task can be 
greatly simplified by transforming to the "group coordinates," defined with respect to the 
basis set composed of vectors which transform according to different irreducible represen- 
tations contained in T, in which the Jacobian is block-diagonal. In practice, it is usually 
impossible to determine whether the isotropy group exhausts the dynamical symmetries 
of the system or the group £ contains some hidden symmetries as well. It is, therefore, 
important to show that a number of restrictions on the set of control parameters can be 
obtained using an arbitrary unitary subgroup CJ of C. 

Decomposing the representation T into a sum of irreducible representations T r of the 
group C with respective dimensionalities d r , we obtain: 

T = piT 1 © p 2 T 2 © • • • © p g T q (28) 

with 

n x =pid 1 + p 2 d 2 -\ Vp q d q , (29) 

where p r denotes the number of equivalent representations T r present in the decomposition 
([28]), and q is the total number of nonequivalent irreducible representations. Since CJ is 
unitary, all irreducible representations T r in ( p8|) can be chosen as unitary [|16j. 

The tangent space T is similarly decomposed into a sum of invariant subspaces Ug such 
that T(g)x E L™, Vx 6 Ug and V# G £': 

T = L l C i © L 2 C , © ■ ■ ■ © L q c ,, (30) 

where 

L £/ = Lj*/ © L/£i © • • • © L £t (*^1) 

and a — 1, • • • ,p r indexes different invariant subspaces, which correspond to the same group 
of equivalent irreducible representations T r . It should be noted that even though the de- 
composition (|30|) is unique, the decomposition (|3lD is not, unless p r = 1. Let us introduce 



S 



a basis in each invariant subspace and denote the basis vectors e™, i = 1, • • • , d r . We 
choose the basis vectors such that they transform according to the irreducible representation 
T\ i.e., 



(32) 



For unitary T r a generalized orthogonality condition between basis vectors e[ Q can be es- 
tablished [T(J as a consequence of 



(ef ■ ef ) 



5 r , s 5 M (ef ■ eD 



(33) 



In addition, for p r > 1 the decomposition (|3~TD can always be performed in such a way that 
(e[ • e[ Q ) = 5 a> /3 (this, however, still leaves some freedom in choosing the invariant subspaces 
L™), so that the complete set of basis vectors {e™}, where r = 1, • • • , q, a = 1, ■ ■ • ,p r , and 
i — 1, • • • , d r is made orthonormal. We, therefore, conclude that the matrix P defined by 

(34) 





-p 1 ' 




' P{ ' 




- (eI i)t - 


p = 




pr 




pr 






. pq . 




pr 




M Pr Y. 


r x = 


pt( 


or, more generally, 


unitary). 





Furthermore, according to the Wigner-Eckart theorem ||16|| , the matrix elements of an 
arbitrary matrix (and the Jacobian A, in particular) invariant with respect to any group 
transformation 



T(g)AT-\g) = A, Vg e £', 
satisfy the following general formula: 



(ef ■ AeV*) 



and the scalar product 



r/3\ 



(35) 



(36) 



(37) 



is independent of the index i 



,d r (but depends on the decomposition (|31|)). As 



a result, on transformation to the group coordinates the Jacobian matrix becomes block 
diagonal: 



^A 1 



A = PAP' 1 



Ai 



(38) 



where each block A r is itself block-diagonal 

A r 



A'' 



(39) 
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and consists of d r identical p r x p r blocks A r with the matrix elements defined by the scalar 
product fl37l). 

If no irreducible representation T r of £ enters the decomposition ( p8| ) more than once, 
i.e., pi = ■ ■ ■ = p q = 1, the structure of the Jacobian matrix is completely resolved: the trans- 
formed Jacobian is diagonal and its spectrum consists of eigenvalues A r = A r , r = 1, • ■ ■ , q 
with multiplicity d r , while the basis vectors e™ become the corresponding eigenvectors (and, 
consequently, define the normal modes of the linearized system). In this case the invariant 
subspaces of the group £ define the eigenspaces of the Jacobian, U = U a . Clearly, the 
spectrum becomes degenerate, if the symmetry is sufficiently high (such that T contains at 
least one irreducible representation T r with dimensionality larger than one). 

Degeneracy should not necessarily be associated with symmetry and might be accidental 
(with respect to the group £,'). For instance, it can happen that A r = A r for some r ^ r', so 
that the multiplicity of the eigenvalue A r is increased respectively to d r + d r r. Accidental de- 
generacies can be alternatively thought of as a consequence of hidden symmetries contained 
in the full symmetry group C of which CJ is a subgroup. However, the degeneracies not 
associated with some physical symmetry are likely to disappear under a typical perturba- 
tion, such as a change of system parameters and, therefore, are most conveniently regarded 
as accidental. Since the full symmetry group C, in general, depends on system parameters 
and cannot be directly deduced from the structural symmetry group Q, it is usually more 
convenient to use its parameter-independent subgroup CJ = Ti^ instead. 



D. Jordan Decomposition 



If the symmetry described by CJ is low, a number of equivalent irreducible representations 
will typically be found in the decomposition (EH), i.e., we will have p r > 1 for certain r. In 
this case the knowledge of the dynamical symmetries alone is not sufficient to completely 
determine the structure of the Jacobian matrix, which is, in general, system-dependent. As 
a result, one has to solve a secular equation 



(A r - \ T J)e r 



0. 



(40) 



for each block A r with p r > 1 in order to find the eigenvectors in the invariant subspace 
U' c , and the respective eigenvalues. Here, unlike the case of quantum mechanics, the Ja- 
cobian matrix does not have to be Hermitian and, therefore, might not be diagonalizable. 
However, A r can always be reduced to the Jordan normal form by finding the coordinate 
transformation Q r such that 



A r = Q r A r (Q r 



A 



rl 



A^ 



(41) 



where p' r < p r is the number of distinct eigenvalues and the Jordan superblock 

'A™ 



A r 



A™ , 

Jra 



(42) 
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corresponding to the eigenvalue X ra consists of j ra Jordan blocks 

X ra 

\ra ^ ^™ 

i 

1 X rt 



(43) 



In the absence of accidental degeneracy all eigenvalues of A r are different, so that p' r = p r 
and j ra = 1 for all a, i.e., A r is diagonal. 

Since each block A r of the transformed Jacobian ( |38"D consists of d r identical blocks A r , 
applying the coordinate transformation defined by the block-diagonal matrix assembled from 
d r blocks Q r , 



Q r 



Q' 



(44) 



reduces A r to the Jordan normal form: 

A r = Q r A r (Q r Y l 



A r 



A' 



(45) 



The Jordan blocks on the diagonal of A r will not, in general, be arranged in superblocks 
with the same eigenvalue. This, however, can be trivially corrected by permuting the rows 
and columns of A r to obtain the matrix 



A r = R r A r( R r\-l 



A r Pr 



(46) 



where R r is the permutation matrix arranging the identical Jordan blocks next to each other, 
and the Jordan superblock corresponding to the eigenvalue X ra has the form 



A r 



A\ c 



Jra 



(47) 



Each block A\ a is, in turn, composed of d r identical Jordan blocks A™, defined by (|43f) : 

-Ara 



A? 



(48) 



Defining the block-diagonal coordinate transformation matrices Q and R 



■Q 1 



R 1 



Q 



R 



(49) 



Q" 



R q 
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we eventually obtain the sequence of coordinate transformations reducing the Jacobian ma- 
trix A to the Jordan normal form: 



A = (RQP)A(RQPy 1 = 
where each block A r , r = 1, ■ ■ • , q is defined by 



A 1 



A q 



(50) 



E. Conditions for Controllability 



Once the Jacobian is reduced to the Jordan normal form, we can turn to the problem 
of reducing the controllability condition to a set of simpler conditions that will give us 
the restrictions on the admissible set of control parameters. Since the controllability is 
a property of the system which does not depend on the choice of the coordinate system, 
condition flT3| ) is invariant with respect to any (nonsingular) coordinate transformation 
and hence is satisfied for the pair (A, B), if and only if it is satisfied for the pair (A,B), 
where B = (RQP)B is the transformed control matrix. Let us partition the transformed 
control matrix B according to the block structure of A: 



B 



B 1 
B q 



B' 



B'P'r 



B\ l 



3ra 



B] 



jbra 



r>ra 
D id r -I 



(51) 



and denote the first row of the matrix B™. Next, define the matrix B ra using the 
relations 



B' 



B[ c 



U id r 



(52) 



In the absence of accidental degeneracy between the eigenvalues that correspond to different 
invariant subspaces L r c , equation P7| ) ensures that there are exactly d r j ra Jordan blocks A[ a 
with the same eigenvalue X ra . If, however, there is such an accidental degeneracy involving 

, L r £,, such that for certain a\, - ■ ■ ,a s 



s different invariant subspaces If^,, 



A 



riai 



A, 



(53) 



the number of Jordan blocks corresponding to the eigenvalue X ra increases to 



3rd 



^2 ^'jr', 



(54) 



The knowledge of the number of Jordan blocks is very useful, since, according to the 



standard result of linear system theory [11[, it ultimately determines the minimal number 
of control parameters. Specifically, it can be shown that the controllability condition for the 
pair of matrices (A, B) is satisfied, if and only if for every r and a (taken equal to r± and 
ai below) 
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rank 



B r 



B' 



Jra 



(55) 



where the indices and aij are chosen according to ([xf). This, in turn, can be achieved, if 
and only if n u > j' ra for every r and a. Hence, in the most general case the minimal number 
n u of independent control parameters should equal the maximal number of Jordan blocks 
with the same eigenvalue \ ra : 



n„ 



max max j ra . 

r=l,---,q a=l,—,p' r 



(56) 



Note that, since the block B ra has d r j ra rows, rank(5' ra ) < d r j ra for every r and a. 
Using this fact, the trivial matrix inequality 



rank 



B' 



< rank(P/ iai ) + • ■ ■ + rank(P/ sas 



(57) 



and equation (|55|) one obtains 

rank(P/ ra ) = d r j ro 



r = !,•••, q, a = l,---,p r - 



(58) 



Furthermore, according to the definition ( p2[ ) of the matrix B ra , for every r and a 
rank(S ra ) > rank(S ra ), so one can write 



rank(-B r ) = rank 



rl n 



B 



> max rank(5 ra ) > d r max j ra . 



a=l,-,p 



a=l, — ,p r 



(59) 



In addition, since Q r and R r are nonsingular coordinate transformations which do not change 
the rank of a matrix, 



rank(iT) = rank(iTQ r P r B) = rank(P r 5), 



(60) 



where we got rid of all system-specific information, which was contained in the matrices Q r 
and R r . 

The symmetry information alone is insufficient to determine the values of either j ra or 
f ra . However, by definition one has p r > j ra > 1 so that f ra > d r . As a consequence, 
we obtain two necessary conditions for controllability. First of all, equation (|56| ) yields the 
lower bound on the minimal number of control parameters 



n u > max d r . 

r=l,-,q 



(61) 



Second, inequality (|59|) combined with equality (|60D imposes a number of restriction on the 
control matrix B, 



rank(P r 5) > d r , r = 1, • • • , q, 



(62) 
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which can be interpreted as the requirement of the mutual independence of control pa- 
rameters. We can therefore, conclude that an arbitrary (unitary) subgroup CJ of the full 
dynamical symmetry group C does not completely define the minimal set of control param- 
eters. It does, however, define a set of necessary conditions required for controllability. In 
general, the knowledge of all dynamical symmetries, both unitary and nonunitary, described 
by the group C is required in order to completely resolve the structure of the Jacobian matrix 
and obtain the necessary and sufficient condition for controllability. 

Nevertheless, even without knowing the full symmetry group C one can obtain the neces- 
sary and sufficient conditions by making a number of assumptions. First, assume that there 
are no accidental degeneracies (it is usually safe to do so if, e.g., CJ is taken to coincide with 
Ti^'. we ensure that all physical symmetries are taken into account, and accidental degen- 
eracies should only appear for certain special values of system parameters). Then j ra = 1, 



j' ra = d r , and B ra = B ra for all r and a, so condition (|56"D is equivalent to 



n u = max d r . (63) 

r=l,---,q 

If, in addition, no irreducible representation T r of £ enters the decomposition ( ^8|) more 
than once, such that p r = 1 for all r, instead of inequality (|62|) one obtains the equality: 



rank(P r P) = d r , r = l,---,q. (64) 



Conditions fl5"2| ) and ( |B"4"D can be simplified even further by defining the projection op- 
erator P r = (P r )tp r onto the invariant subspace U a C T. This operator can be obtained 
directly from the matrix representation T for most symmetry groups of interest. For finite 
discrete groups it is given by 

P r = ^J2x r (g)Ag), (65) 
9 gee 

where n g is the number of elements of the group £ and x r {g) ls t ne character of the group 
element g in the representation T r . Similarly, for compact continuous groups we have 

P r = d r [ x r (9)T(g)d^g), (66) 

where dfi(g) is the group measure ||16| . Observing that rank( (P r y P r B) = ia.nk(P r B), we 
can use the projection operators to rewrite the condition (]64]) in an equivalent form 



rank(P r 5) = d r , r = 1, • • • , q. (67) 

Summing up, we conclude that with the two assumptions made above the system is 
controllable, if and only if the two conditions are met. The first one requires the number 
n u of control parameters to be greater or equal to the dimensionality d r of the largest 
irreducible representation T r present in the decomposition of the matrix representation T of 
the subgroup £ C C in the tangent space T. The second one requires the control parameters 
to be independent: the columns h m of the control matrix B have to be chosen such that 
d r of the projections P r h m , m = 1, • - - ,n u are linearly independent (and, therefore, span 
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the eigenspace U = L r c ,) for every r = 1, ■ ■ ■ , q. The last requirement imposes a number of 
restrictions on the admissible form of the linear response of the system to perturbations of 
control parameters. 

A number of comments are in order. First of all, as we have just seen, the number of 
control parameters is determined by the number of Jordan blocks with the same eigenvalue, 
not the multiplicity of that eigenvalue. It becomes intuitively clear why this is so, if one 
compares the action of different Jacobians already reduced to the Jordan form. For instance, 
the Jacobian 



A 1 



A 



A 



(68) 



generates the set of three linearly dependent vectors h° = b, h 1 = Ab, h 2 = A 2 b (compare to 
(|TT|)), that span a one-dimensional subspace of M 3 for an arbitrary choice of b. As a result, 
three control parameters and a control matrix with three linearly independent columns, 
B = [h% b 2 b 3 ], are necessary to control the system. On the contrary, the Jacobian 



A 

1 A 
1 A 



(69) 



generates a linearly independent set of basis vectors that spans M 3 , requiring just one control 
parameter and a control matrix with a single column B = b. 

Second, symmetry does not always make the Jacobian degenerate, and the nondegenerate 
case can be handled in the same way as the one with no symmetries. Neither does the 
degeneracy by itself imply that multi-parameter control is required: even if the eigenvalue 
X r ' a ' is degenerate, but f ra — d r — 1 for every r and a (the degeneracy is accidental and 
limited to a single invariant subspace L r c ,), one control parameter is sufficient to ensure 
the controllability. In both cases, however, the dynamical symmetry should be rather low. 
Specifically, the decomposition fl28|) of the matrix representation T should not contain any 
multi-dimensional irreducible representations. 

Finally, the conditions on the set of control parameters that were obtained above are im- 
posed by the controllability condition and guarantee that control can be achieved. However, 
in general, only the weaker stabilizability condition has to be satisfied which, according to 
section fl A] , requires that every unstable normal mode of the system is controllable, so that, 
only r and a such that | X ra \ < 1 have to be considered in the conditions (|55| ) and (|56D . As a 
consequence, it might be possible to stabilize highly symmetric states of compact extended 
systems with strong spatial correlations using a single control parameter — if only a small 
number of modes is excited, there is a chance that all unstable modes will correspond to 
one-dimensional irreducible representations T r . In strongly chaotic systems a large num- 
ber of modes will be unstable and many of them will inevitably be degenerate, calling for 
multi-parameter control. Similar considerations apply to weakly chaotic systems with large 
spatial extent. 
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III. TIME-PERIODIC STATES 



The results obtained above for the time-invariant case can be generalized for the time- 
varying and, in particular, time-periodic case, but first we have to define the notions of 
controllability and dynamical symmetry in the context of time-varying trajectories. Indeed, 
in the time- varying case the Jacobian A 1 and the control matrix B l in the linearized evolution 
equation (|]) are time- dependent and, as a consequence, neither the definition of controlla- 
bility given in section |11 A| nor the condition ( |T3| ) holds. Besides, it is not at all clear that 
the symmetry of the target trajectory, and hence the dynamical symmetry group £ can be 
uniquely and consistently defined. 

We will see that all these notions generalize in a rather straightforward way, so that 
the same formalism as we used in the previous sections applies here as well. To begin 
with, we define the controllability of a general time-varying linear system. Expanding the 
definition given for time-invariant target states, we call the ^^.-dimensional linear system 
(0) or the sequences of matrices {A f , B 1 } controllable if, for any initial state Ax* 1 = Ax.;, 
times tf — ti > n x , and final state Axj, there exists a sequence of control perturbations 
Au*% • • • , Au*-^ 1 such that the solution of equation (0) satisfies Ax^ = Axj. 

The controllability condition can be restated in terms of the matrices A 1 and B l con- 
ducting the analysis similar to that of section |11 A| . Applying the map (Q) n x times yields 

Tlx— 1 

Ax t+n. = j£+n*-l Ax t + J*+^:^* +fc Au* +fc , (70) 

k=0 

where we have introduced a shorthand notation 

J* = A t A t ~ l ■ ■ ■ A l - k+1 (71) 

for the product of k consecutive Jacobians. Arguments identical to those used to derive the 
controllability condition ([13]) from equation (|) allow us to conclude that for time-varying 
states the controllability condition can again be written in the matrix form: 

rank(C t ) = n x , Vt, (72) 

where the controllability matrix (|14|) is now replaced with the sequence of matrices 

Ct=[B* J^B*- 1 4B 1 - 2 ■■■ J^B 1 -^ 1 }. (73) 

Next we have to define the dynamic symmetry group C. Suppose the target trajectory 
x^x 2 , • • • ,x T has period r, and the symmetry of the point x* on the target trajectory is 
described by the group Ti^t C Q. We can then write 

S(x* +1 ) = «7(F(x*, u)) = F(0(x*), u) = F(x*, u) = x' +1 (74) 

for every g G Ti^t. Consequently, 

T^x 1 Q 7~(-x 2 Q " ' Q 7~Lyj Q T^* 1 , (75) 
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which means that the symmetry properties of all the points on the target trajectory are the 
same and the isotropy symmetry group of the trajectory TC^ can be uniquely defined using 
an arbitrary point x*, TL^ = H^t. 

Using the arguments that lead to equation (|20|) we obtain for an arbitrary g e Ti^: 

x <+1 + #(A'Ax) = #(x <+1 ) + #(A'Ax) = #(x m + A 1 Ax) 

= s(F(x*, u) + A* Ax.) = g(F(9* + Ax, u)) 

= F(<?(x« + Ax), u) = F((/(x') + p(Ax), u) 

= F(x* + g(Ax), u) = F(x', u) + A'#(Ax) 

= x m +^(Ax). (76) 

This, in turn, means that the symmetry group C t of the Jacobian A 1 satisfies 

WxC£ t) t = l,---,r. (77) 

Again, typically, we expect C t = H^, so that C too would be unique for any given periodic 
trajectory as would the matrix representation T, such that 

T(g)A t = A t T(g), \/g G C. (78) 

It is, therefore, enough to know the symmetry properties of an arbitrary point of the periodic 
trajectory in order to establish the requirements on the control scheme similarly to the time- 
invariant case. If C t is not unique, we can still use the commutation relation ([7|) for the 
subgroup £' = 7ix to obtain a lower bound on the minimal number of control parameters. 

Finally, we note that although it is possible to obtain certain results for time-varying 
control matrices B t , we assume, as is often the case in real systems, that B l is constant and 
drop the time index. As we will discover below, in the time-periodic case the restrictions 
imposed by symmetry on the structure of the matrix B can typically be determined without 
the detailed knowledge of the Jacobian matrices, but based on the symmetry properties 
alone, similarly to the time-invariant case. Indeed, let us construct the representation T of 
the group CJ in the tangent space T and decompose it into the sum of irreducible repre- 
sentations. This again defines a set of invariant subspaces L™ and a set of basis vectors 
{e[ Q }, which we use to construct the coordinate transformation matrix P according to the 
definition (||). 

Since the rank of the matrix (|7^) does not change under a coordinate transformation, 
the controllability condition ( |72| ) is equivalent to the condition 

rank(C t ) = n x , t = 1, • • • , r (79) 

where 

C t = [B J\B ••• Jl^B], (80) 

B = PB and J\ = PJ^(P) -1 . The products J l k have the same symmetry properties as the 
Jacobian matrices A 1 for arbitrary k and t, and, therefore, both the matrices A 1 and the 
products J l k block-diagonalize in exactly the same way: 

A t = PA t (P)~ 1 = ••• , (81) 

fit* 
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and 



Jl = PJtiP)- 1 



J, 



(82) 



J k -I 



Similarly to the time-invariant case, the blocks A t,r and j]f are themselves block- diagonal 



A 



t,r 



1 fc 



J 



1 fc 



(83) 



and consist of d r identical p r x p r blocks A*' r and r£ r , respectively, whose matrix elements 
are defined by the scalar products 



— (o ra A t a r l 3 \ 

afi = \£i -A e i ), 

J- k )afi= {^i ■ ). 



(A*' r ) 



(84) 



Using the definition ( |7I| ) one can check that for any t, k and r the matrix r fc ' r can be 
represented as the product 



-p<,r- p^t,r j^t—l,r jyt—k+l,r 



(85) 



Let us partition the transformed control matrix B into blocks B r = P r B and define the 
reduced controllability matrices 



c; r>" -i' ; r n r ■■■ ./;;" ; IV . (86) 

Using relation (|2^) and the fact that the matrix C[ has d r p r rows one can write 



rank(C t ) = rank 



C 1 



C q 



< rank(C t 1 ) + • • • rank(C^) < dxp\ ■ ■ ■ d q p q 



(87) 



to obtain as a consequence of ( [T5[ ) the set of reduced controllability conditions 

rank(C[) = d r p r , r = 1, • • • , q. (88) 

The blocks B r of the matrix ( j36"l) can become linearly dependent for certain r, d r 
and p r . Indeed, it is trivial to see that for a sequence of n arbitrary p x p matrices i?,, it is 
always possible to find a set of coefficients fi , fix, • • • , // n such that 



fl I + fllRl + ■ ■ ■ + fi n Rn — 0, 



(89) 



as long as n > p 2 . Equally easy to establish is the fact that, if the matrices Ri are not 
arbitrary, but satisfy the condition 



R l = W 1 W 2 ---W l 



(90) 
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where Wi is a sequence of arbitrary p x p matrices, such that Wi +T = Wi, equation ( p9|) can 
always be satisfied for n > min(p 2 ,pr). The p r x p r matrices fy, • • • , f^_i form precisely 
the sequence satisfying the condition fl9~0|). Besides, if the condition fl89|) is satisfied for 
Ri = r*' r , it is satisfied for the sequence Ri = j\' 7 as well. As a result, 

rank(C[) = rank[5 r ?{ r W ■■■ T n \B r \ (91) 

for n = min(p^,p r r) and arbitrary B r . Therefore, in order for the conditions flS~8D, and hence 
(|72]), to be satisfied, one should have 

rank(P r E) > ceil ^max ^j^j , r = 1, • • ■ , q, (92) 

where ceil(x) denotes the smallest integer number n such that n > x. The necessary condi- 
tions on the control matrix B, defined by (^) are the generalization of the time-invariant 
result (|62|). Instead of (61) one respectively obtains the restriction on the minimal number 



of independent control parameters required to satisfy the controllability condition (^) for 
a periodic target trajectory: 

n u > ceil ( max max (—,— ]) • (93) 

\r=l,-,« \p r T J J 

It is interesting to note that a periodic trajectory can be made controllable using the number 
of control parameters n u that could be smaller than the number required for a steady state 
with the same symmetry. 

Three special cases deserve separate consideration. First of all, suppose that the Jacobian 
matrices A 1 commute with each other, so they can be simultaneously diagonalized. In this 
case the condition ([89]) can be satisfied by an appropriate choice of coefficients fii, ■ ■ ■ , n n for 
n > p r , so the necessary conditions fl9"3] ) and ( §2 ) will reduce to (pip and (§^), respectively, 



and n u will no longer depend on the period r of the target trajectory. 

Next, suppose there are no accidental degeneracies between the eigenvalues of the Ja- 
cobians A 1 and their products J l kl and no irreducible representation of £ appears in the 
decomposition ( ]2lf ) more than once (so that Jacobian matrices can again be simultaneously 
diagonalized). Now, however, identically to the time- invariant case one obtains the necessary 
and sufficient conditions (|63| ) and (|64D instead of the necessary conditions (|6l| ) and (|62|). 

Finally, although we used the fact that the trajectory is periodic to derive the above 
results, this requirement could be lifted, provided the symmetry of all points on the target 
trajectory is the same, and, therefore, the condition (|78| ) is satisfied. A nonperiodic trajec- 



tory could then be treated as a periodic one, with period r = oo, and the condition (p9[) 
will be satisfied by an appropriate choice of coefficients fix, ■ • • , \i n for n > pi. As a result, 
instead of the restriction (|93| ) one will obtain 

n u > ceil ( max — J . (94) 



r=l, 



,q p T 
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IV. CONTINUOUS-TIME SYSTEMS 



Most of the results obtained in the previous sections can be directly and naturally gen- 
eralized to continuous-time systems. This is a rather valuable asset of the developed theory, 
since continuous-time control is, in general, a much more flexible and powerful technique 
than discrete-time control. In the presence of a decent continuous-time mathematical model 
(H), continuous-time control can often achieve far superior results. It is, however, a much 
more complicated technique as well. For simplicity we only discuss the control of time- 
invariant target states. Linearizing the evolution equation fl5|) around the steady target 
state s, one obtains 

As(t) = AAs(t) + BAu(t), (95) 
where similarly to the discrete-time case we define the Jacobian 

A = D s $(s,u) (96) 

and the control matrix 



fi = D u $(s,u). (97) 

The symmetries of the nonlinear evolution equation fl5|), the target state s, and the lin- 
earization fl9~5]) are determined identically to the discrete-time case using the relations (|TSD, 
(|19l), and (|2T|), yielding the symmetry groups Q, Hs, and C, respectively. The definitions of 
the notions of stabilizability and controllability in the continuous-time case are completely 
analogous to the ones given in section [11 A| for the discrete-time case. 

The dynamical system described by equation (|9~5|) or the pair (A, B) is said to be control- 



lable if, for any initial state As(U) = Asj, times tf — ti > and final state Asf, there exists 
a (piecewise continuous) control perturbation Au(£) such that the solution of equation fl9~5"|) 
satisfies As(fy) = As/. Otherwise the system or the pair (A,B) is called uncontrollable. 

Similarly, the dynamical system or the pair (A, B) is said to be stabilizable, if there exists 
a state feedback Au(£) = — KAs(t) making the system stable, such that all eigenvalues of 
the matrix A' = A — BK have a negative real part, Re(A' fc ) < 0, Vk. Otherwise the system 
or the pair (A, B) is called unstabilizable. 

The controllability of the pair (A, B) again ensures that all eigenvalues of A' can be 
chosen appropriately, so that any controllable continuous-time system is stabilizable as well. 
The controllability of a continuous-time system is also established using the same criterion 
(|13|) used to test for the controllability in the discrete-time case. As a result, the conditions 
imposed on the control matrix B by the controllability condition in the presence of symmetry 
are exactly the same as those obtained for discrete-time systems. 



V. SYMMETRY VIOLATION 



In reality symmetries of physical systems displaying dynamical instabilities are almost 
never exact. Indeed, the cylinders in a Taylor- Couette experiment are never perfectly cir- 
cular, the temperature inside a chemical reactor is never absolutely uniform, neither are 
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the rotor blades of a turbocompressor exactly identical. The above analysis, on the other 
hand, has been conducted in the assumption of exact symmetry. Therefore, it is essential 
to understand how the obtained results change, if the symmetry is not exact or, in other 
words, what the effect of a weak symmetry violation is. Such an analysis is also crucial 
in the vicinity of points in the parameter space where symmetry increasing bifurcations or 
accidental degeneracies occur. 

For simplicity let us again consider the time-invariant case. The Jacobian A of a weakly 
perturbed symmetric system takes the form 

A = A + eA 1 , (98) 

where e denotes the magnitude of the perturbation and the unperturbed Jacobian A is 
exactly symmetric with respect to all transformations g of the group C. For the group 
representation T we thus have 

T(g)A - A T(g) = 0, \Jg e C. (99) 

In general, the perturbation eA\ will not be symmetric with respect to any element of the 
group £, except the identity transformation e: 

T{g)A l -A 1 T{g) 1 LQ, Wg G £ \ {e}. (100) 

Therefore, since 

T(g)A - AT(g) = e(T(g)A 1 - A 1 T(g)), (101) 

the perturbation (|9~8"D completely destroys the symmetry of the linearized evolution equation 
([]) for any e ^ 0. As a result, the perturbed system can be made controllable using a single 
control parameter, irrespectively of the properties of the original symmetry group C. For 
instance, calculating the controllability matrix of the perturbed system with n u = 1 and 
B = b one obtains 

C = C + ed + o(e 2 ), (102) 

where we defined 

C = [b A b ■■• (i4 ) n »- 1 b], 

Ci = [0 A x h ■■■ ((Ao) n °- 2 A l + --- + A 1 (A ) n *- 2 )b}. (103) 

Co is clearly the controllability matrix of the unperturbed system with full symmetry, which 
does not have a full rank, if the decomposition ( f2l|) contains at least one irreducible represen- 
tation T r with the dimensionality d r > 1. Indeed, in the absence of accidental degeneracies 
that would mean 

<i 

n = rank(Co) < ^^£V < n x . (104) 

The controllability matrix C of the perturbed system, on the other hand, has full rank for 
any e / because the symmetry is completely destroyed by the perturbation. Therefore, 
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the perturbed linear system becomes controllable even though the unperturbed system is 
not, for arbitrarily small perturbations. 

The controllability ensures that for any initial and final states of the linear system (||) 
the control can be found mapping the initial state to the final state in n x iterations. Using 
(0) one obtains explicitly 



AU* = 



Am* 



(C)- L (Ax t+n * - (A) n *Ax*). 



(105) 



Formally, if the system is controllable, the controllability matrix is invertible, and the so- 
lution ( |105| ) is well defined for any Ax* and Ax* +n *. However, when the matrix C is close 
to being singular its inverse is not well defined. It is convenient to use the singular value 
decomposition of the controllability matrix 



where Q = [q a q 2 
matrices, and 



C = QZR\ 
q nx } and R = [ r 1 r 2 



(106) 



r nx ] are some orthogonal n x x n x 



(107) 



The singular values are ordered such that 0i(e) > ^(e) > • • ■ > (Tn x { e ) ^ OT Additionally, 

(108) 



equation (|104|) requires 



lim cTi(e) = 0, % = n Q + 1, ■ ■ ■ , n x . 

e^0 



In terms of the matrices Q, E, and R we can write the inverse of C as 

(C)- 1 = i?(E)- 1 Qt = ^ (T ri (e)riq t 



(109) 



i=l 



and, therefore, for small e equation (|105|) gives 



AU* ~ ^ ,Axt+W ')-^-^ W ' Axt ) r , 



«=no+l 

As a consequence, we obtain 



lim|AU*| 



oo. 



(110) 



111) 



This relation means that at least one control perturbation of the feedback sequence 
Au*, • • • , Au t+rix ~ l diverges as the symmetry breaking perturbation eA\ of the Jacobian 
vanishes. Since no specific relation between the initial and the final state of the system was 
implied, the obtained result is general, and does not depend on the control method used to 
calculate the feedback. 
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In fact, a more general statement holds. Suppose the symmetry is violated only par- 
tially, such that the perturbed Jacobian (|98|) remains exactly symmetric with respect to a 
subgroup CJ of the full symmetry group C. Denote n u and n' u the minimal number of control 
parameters required (assuming exact symmetry) by the groups C and C, respectively. Then 
it can be shown that, similarly to the single-parameter case, at least one control perturbation 
of the feedback sequence Au*, • • ■ , Au t+nx ^ 1 diverges as the symmetry breaking perturba- 
tion eAi of the Jacobian vanishes whenever n' u < n u < n u . The same result is obtained if 
the independent with respect to the group CJ control parameters become dependent with 
respect to the group C, as indicated by the violation of the general independence condition 
(|55|). The time-periodic generalization is also straightforward. We will call this situation 
parametric deficiency. 

In other words, although it might be possible to control a linear system with approximate 
symmetry using a number of control parameters which is smaller than that required in the 
assumption of exact symmetry, the stabilization requires feedback of very large magnitude. 
Such systems are called weakly controllable in the language of control theory. However, the 
linear system is only an abstraction. The linear approximation (0) of the evolution equation 
(H) is only valid for small perturbations Au* of the control parameters and small deviations 
Ax* from the target trajectory. Besides, additional restrictions on the magnitude of the 
feedback are usually imposed by practical limitations, size and energy constraints, etc., at the 
implementation stage. One can, therefore, conclude that, since the feedback scales linearly 
with the deviation from the target trajectory, a nonlinear system with parametric deficiency 
can be stabilized using linear control only in an asymptotically contracting neighborhood of 
the target trajectory. 

Finally, consider the vicinity of the point u in the parameter space M ra " at which an 
accidental degeneracy occurs, such that the dynamical symmetry is described by the group 
£ for u 7^ u and is increased to C (of which CJ is a subgroup) for u = u . In this 
case C can be considered approximate symmetry in the vicinity of uo, and the distance to 
that point determines how strongly (or weakly) the symmetry C is violated. Suppose the 
control scheme is such that there is a parametric deficiency. Then the system will remain 
controllable for However, the strength of feedback required to control the system 

will diverge as u approaches u , at which point the system will become uncontrollable. 

VI. SYSTEM IDENTIFICATION 

Finally, let us consider the problem of phase space reconstruction (sometimes also called 
system identification) in the context of symmetric systems. The need to use phase space 
reconstruction stems from the fact that few experimental systems exist whose dynamical 
equations are known (let alone with precision required by the control algorithms). Even the 
complete information about the state of these systems is rarely available, so one typically 
has to contend with having a measurement of a single scalar output (sometimes called an 
observable) for the description of the dynamics. This output, which we call y(t) is, in general, 
a function of the (unknown) internal state s(t) of the system described by the (also unknown) 
continuous-time evolution equation (^). In other words, 

y(t) = G(s(t)). (112) 
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It turns out that it is possible to reconstruct both the internal state of the system and 
its dynamics based on the time series measurement of the output y(t) using the procedure 
originally proposed by Packard et al. fIT|j . The easiest way to recreate the state of the 
system from the scalar output signal is to use time delays. Let us choose different delay 
times Ti, T 2 , • • ■ , T nx and construct an n^-dimensional delay coordinate vector 

v(t+n) 

y(t + T 2 ) 



x(f) 



(113) 



ly(t+T nx 

Takens showed \TB\ that for a scalar output ( |112j ) and conveniently chosen delay times T i; if 
the dimension n x of the embedding space is such that n x > 2n^ + l, where n 1 ^ is the Hausdorff 
dimension of the attractor A, the map P : s(t) — > x(t) generically provides a global one-to- 
one representation of the attractor and, hence, the system state. Finally, one recovers the 
map (|I|) by defining x fe = x(t fc ), where t , t\, t 2 , ■ • • denote the times of consecutive crossings 
of a conveniently chosen Poincare surface of section in the embedding space M. nx by the 
reconstructed trajectory x(i). 

However, although the Takens' embedding theorem holds for a typical system without 
symmetries, it is not satisfied for most of the symmetric systems due to the symmetry- 
related degeneracy of the evolution operators. In other words, most symmetric systems are 
nongeneric in the sense of Takens. As a result of this nongenericity the state of the system 
becomes impossible to reconstruct using a single scalar output, no matter how large the 
dimensionality n x of the embedding space is, locally in the vicinity of highly symmetric 
periodic trajectories. The attractor of the system remains folded at certain points and along 
certain curves in the phase space, which prevents the global reconstruction as well. Using 
the language of control theory we will say that such systems are unobservable locally as well 
as globally. 

The question of symmetry-caused nongenericity in the framework of phase space recon- 
struction of a general symmetric system was first considered by King and Stewart [Tj| , who 
determined that the reason for the failure of the embedding theorem is the violation of one 
of Takens' assumptions that the flow defined by equation (|5|) has simple eigenvalues for 
low-period periodic trajectories. As we have seen above, symmetric systems typically (but 
not always) have degenerate eigenvalues (due to the fact that most nontrivial irreducible 
representations are multi-dimensional) and, as a consequence, are nongeneric. King and 
Stewart went on to formulate and prove a generalization of the Takens' embedding theorem, 
which required the output to be a vector, not a scalar, function of the actual state of the 
system s(t): 



y(t) = G(s(t)), 



(114) 



mapping the phase space Q of the original system onto an n^-dimensional Euclidean space. 
The state of the system can then be similarly represented by a delay coordinate vector 

y(f + ?l) 
y(* + t 2 ) 



x(t) 



y(t + T n 



;ii5) 
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where now the dimensionality of the embedding space is n x = n y n e . The question we have 
to answer is what conditions should the function G satisfy in order to allow a local (or 
global) one-to-one embedding. Since the exact form of the evolution equations is assumed 
unknown, in order to find the answer one can only exploit the symmetry properties of the 
system, which are often easy to establish based on the underlying symmetries of the physical 
space. Fortunately, the symmetry provides most of the necessary information. 

Since we are interested in the issue of phase space reconstruction only as far as it applies 
to the problem of linear control, we will assume a local character for the observability 
property, unless explicitly stated otherwise. According to the analysis conducted in [T3|], local 
embedding in the vicinity of the periodic trajectory s(t) requires M. ny to contain at least one 
copy of every invariant subspace L™ generated by the (unitary) irreducible representation 
T r of the respective isotropy symmetry group CJ = He- This would lead one to assume that 
the minimal dimension n y of the output signal should be determined by the dimension of 
the largest irreducible representation T r . 

This assumption can be trivially verified using the formalism developed above for the 
control problem in the presence of symmetry. Indeed, let us again consider a time-invariant 



target state s. Linearizing the output ( 114 ) in the vicinity of this target state and denoting 



the displacement Ay(t) = G(s(£)) — G(s) one obtains: 

Ay(t) = CAs(t), (116) 
where the constant matrix C is defined thus: 



C = D s G(s). (117) 

The dynamical system defined by equations ( pop and ( |116| ) or the pair (A, C) is said to be 
observable if, for any times tf — ti > 0, the initial state As(tj) = Asj can be determined from 
the measurement of control perturbation Au(i) and output Ay(i) in the interval t G 
Otherwise, the system or the pair (A, C) is said to be unobservable. 

It can be shown that the notion of observability is dual to the notion of controllability. 
The crucial benefit of this duality is the fact that the observability condition for the pair 
(A,C) is equivalent [JTTJ] to the controllability condition for the pair (A\&). Since the 
commutation relation (p3j) directly implies that 

T( g y^ = A^T(g)\ VgEC\ (118) 

the symmetry properties of the matrices A and A' are essentially identical (as are the 
structures of their spectra, Jordan normal forms, etc.). As a result, all restriction imposed 
on the control matrix B by the controllability condition in the presence of symmetry should 
be satisfied for the matrix as well. 

For instance, in case there are no accidental degeneracies and the representation T con- 
tains at most one copy of each irreducible representation of the group one can claim 
that in order to reconstruct the dynamics in the vicinity of the time-invariant symmetric 
target state s the number n y of measured scalar output signals y.i(t) should be the same as 
the minimal number n u of independent control parameters, i.e., n y = n u . Furthermore, the 
outputs have to be independent, so that d r of the projections P r Cj, i = 1, • • • , n y are linearly 
independent for every r, where P r is the projection operator defined by (|65D and (|66"D, and 
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c] is the ith row of the matrix C, which imposes a number of restrictions on the allowed form 
of the function G. Local observability can be defined in the same way for discrete-time sys- 
tems. Careful consideration shows that symmetry produces similar constraints independent 
of the particular description. 

In addition to the notion of controllability of individual eigenvectors it is often also 
convenient to define the notion of their observability. We will say that the eigenvector e of 
the Jacobian A is observable, if there exists i, 1 < % < n y , such that (e-Cj) 7^ 0. Respectively, 
an eigenvector that is orthogonal to every row of the matrix C is called unobservable. Clearly, 
the observability of the linearized system is equivalent to the observability of each and every 
eigenvector of the Jacobian matrix. 

In the conclusion of this section we make a few comments regarding the problem of 
global phase space reconstruction. Often it is important to know how the symmetry of the 
continuous-time experimental system transpires in the structure of the discrete-time map 
([[]) obtained as a result of the time delay embedding produced by a general output signal 
( ITTi ). King and Stewart |3| recognized that it is as important to preserve the symmetry 



of the attractor as it is to preserve its topology during the reconstruction. According to 
( |114j) , using an arbitrary vector output y(t) to generate the delay coordinate representation 
of the system state corresponds to picking a function G which, in general, distorts the 
symmetry. In order to preserve the symmetry of the original attractor the function G has 
to be ^-equivariant: 

G(<7(s)) = <?(G(s)), \/g e G, Vs e Q, (119) 

where Q is the structural symmetry group of the system ([5]) which will, in general, act 
differently in the phase space Q and Euclidean space M. nv . In addition, the dimensionality 
n y of the Euclidean space has to be chosen high enough to avoid local folding (obviously, n y 
should be no smaller than the number n y evaluated for the steady or periodic trajectory with 
the highest isotropy symmetry). Finally, a global one-to-one embedding can be achieved by 
choosing n e > 2n^ + 1 to preserve the topology of the attractor. The map ([]]) constructed 
using this embedding will preserve all dynamical symmetries of the original system. However, 
the structural symmetry of the differential equation @ and the map ([!]) will, in general, be 
different. 



VII. EXAMPLES 



A. Coupled Map Lattice 



Now that the formal theory is constructed, we can illustrate it by applying to a few 
simple symmetric systems. First, consider the deterministic coupled map lattice (CML) 
with nearest neighbor diffusive coupling j0|, which is described by the following evolution 
equation 



x 



t+i 



ef(x t t _ l ,a) + (l-2e)f(x t i ,a)+ef(x t i 

+1; a ) 



(120) 



Coupled map lattices of this type are often used as models for such spatiotemporally chaotic 
phenomena as surface growth, population dynamics, and turbulence. Here the index i = 
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1,2, ■ ■ ■ ,n x labels the lattice sites, and the periodic boundary condition is imposed. The 
choice of the map function f(x, a) is usually motivated by the local dynamics of the physical 
system under consideration. In principle f(x, a) can be chosen as an arbitrary (nonlinear) 
function with parameter a, which typically represents the process of generation and growth 
of local fluctuations, while diffusive coupling typically plays the opposite role of dissipating 
these fluctuations. Therefore, the parameters a and e specify the degree of instability and 
the strength of dissipation in the system, respectively. For the purpose of control, however, 
details of the local map are not important. The only aspect of the control problem affected 
by any particular choice is the set of existing unstable periodic trajectories. 

The analysis of the controllability condition conducted in sections and |TTT| shows that, 
if the system is symmetric, certain symmetry-imposed restrictions on the choice of control 
parameters should be satisfied in order to achieve control. In fact, the coupled map lattice 
is by construction highly symmetric. The symmetry is that of the spatial lattice: the 
evolution equation ( |12(J| ) is invariant with respect to translations by an integer number of 
lattice sites (periodic boundary condition makes the group finite) and reflections about any 
site (or midplane between any adjacent sites), which map the lattice back onto itself without 
destroying the adjacency relationship. The corresponding point group C nxV (we assume n x 
- even) has a total of n x /2 + 3 nonequivalent irreducible representations. The first four are 
one-dimensional, d\ = d 2 = d$ = rf 4 = 1, while the rest n x /2 — l are two-dimensional, d r = 2, 
r > 5. In comparison, breaking the reflection symmetry reduces the group to C nx , which 
only has one-dimensional irreducible representations. 

The dynamical symmetry group can be trivially obtained using the observation that 
the Jacobian matrix in the linearization (fj) constructed for the CML ( |120| ) can always be 
represented as a product of two matrices, A 1 = MN f , where 

M,l = (1 - 2e)5 M - + e(S itj ^ + 5 W ) (121) 

describes diffusive coupling, and 

,V/ ; <),./(./,. a )<i,. s (122) 

defines the strength of local instability, with 5ij±\ extended to comply with periodic bound- 
ary condition. This partition of the Jacobian explicitly shows how the symmetry group £ 
depends on the symmetry properties of the nonlinear evolution equation ( |120| ) and those 
of the controlled state x*. The matrix M has all the symmetries imposed by the chosen 
inter-site couplings of the nonlinear model: 

T(g)M = MT(g), Wg G Q, (123) 

while the matrix N f has all the symmetries of the target state x*: 

T(g)N t = N t T(g), V fl GK„ (124) 

where, T denotes the matrix representation of £, and Ti^ is defined as a subgroup of Q. 
Since the Jacobian A 1 commutes with all matrices that commute with both M and N f , we 
conclude that generically £ = Ti^ C Q, in agreement with the general result (|77|). 

Since the symmetry analysis conducted above is valid for every subgroup £' of the dy- 
namical symmetry group, we take £' = £. Constructing the n a .-dimensional representation 
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T of C and decomposing it into the sum of the irreducible representations of C nxV we eas- 
ily determine the restrictions imposed by the symmetry on the minimal number of control 
parameters n u and the structure of the control matrix B. For instance, a zigzag state gives 
C = C nv with n = n x /2 and, according to (|63|) , n u = 2; a non-reflection-invariant state with 
spatial period s corresponds to C = C n with n = n x /s and n u = 1, etc. 

Let us consider the uniform target state, which has the highest symmetry possible, 
C = C UxV , in more detail. The decomposition (ESf) gives 



rpn x /2-\ 



(125) 



and the corresponding basis of normal modes which transform according to these irreducible 
representations is given by the eigenvectors of the operators of translation and reflection, 
i.e., Fourier modes g l : 



cos( jki + 4>i). 



(126) 



Here <pi are arbitrary phase shifts, and ki are the wavevectors defined thus: k\ — 0, ki — 
ki + i = ni/n x for % = 2, 4, 6, ■ • •, and, for n x - even, k nx = n. Fourier modes with the same 
wavevectors k define invariant subspaces L k C T = M. nx . The subspaces L k with < k < n 
correspond to the representations T r with r > 5, L° corresponds to T 1 , and L n to T 4 . 
Since the two-dimensional irreducible representations are present in the decomposition ( |125| ), 
n u = 2. Therefore, in order to control an unstable uniform steady state of the coupled map 
lattice we need at least two control parameters. This is the reflection of symmetric coupling 
in the model (|120|) . Note that, since every two-dimensional irreducible representation occurs 
in the decomposition ( |125| ) once, p 5 = • • • — p nx /2+3 = 1, according to the results of section 
ID, the minimal number of control parameters remains the same for a spatially uniform 
target trajectory of arbitrary time period r. 

On the other hand, since for any length n x of the lattice the group Q = C nxV only 
has one- and two-dimensional irreducible representations and £ is a subgroup of Q, it is 
sufficient to have just two control parameters to make the dynamics of the coupled map 
lattice controllable in the vicinity of a target state with arbitrary symmetry properties 
and temporal period. Choosing the minimal number of control parameters, n u = 2, we can 
determine the conditions making them independent with respect to a particular target state: 
the linear response of the CML to perturbation of the two parameters, given by the columns 
of the control matrix B — [h% b 2 ], has to satisfy conditions ( p7|) and (|6^). 

Failure to satisfy the necessary condition (^7]) rules out the possibility of using global 
parameters, such as the coupling e or parameter a of the local map f(x, a) for control of 
symmetric steady states. Taking u = (a, e), so that 



bx = d a F(5c, 0, u) = M 



d a f{xx,a) 
daf(x nx ,a). 



(127) 



and 



b 2 = <9 £ F(x, 0, u) = (e) (M - /) 



f(xx,a) 
f(x nx} a) 



128) 
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we observe that condition (|27|) is only satisfied, if the group £ is trivial, £ = {e}. This 
result holds for time-periodic symmetric target states as well. 

Alternatively, one can make the system controllable by directly perturbing the system 
at the sites i\ and i<i . The positions of the "control" sites cannot be chosen arbitrary, 
again due to symmetry However, if the target state is spatially uniform, it is trivial to 
show that choosing, e.g., i\ = I and i% = I + 1 satisfies the controllability condition for an 
arbitrary length of the lattice n x . The control matrix corresponding to this choice of control 
parameters can be written in the form = dj iS^i + ^2^,2+1- 

Such localized control also has its downside. In the weak coupling limit, e — > 0, the cou- 
pled map lattice with local feedback becomes a weakly controllable system. The symmetry 
of the lattice of uncoupled maps is described by the permutation group Q = S nx , while the 
linearization about a uniform target state increases the symmetry to £ = GL(n x ): since 
the respective Jacobian is a multiple of the unit matrix, = d x f(x,a)5ij, the linearized 
system is symmetric with respect to all (complex) nonsingular coordinate transformations. 
When coupling is restored, e > 0, the symmetry of both the nonlinear evolution equation 



( |120| ) and its linearization (§) reduces to Q' — £' — C nxV . 

The matrix representation T of the group GL(n x ) in T is already irreducible. Conse- 
quently, n u = n x independent control parameters are required to control the steady uniform 
state of the uncoupled lattice. This result is rather intuitive. Obviously, one can no longer 
control the system applying control perturbations at just two lattice sites, i\ and z 2 . Since 
the control perturbation does not propagate to adjacent sites of the lattice, feedback has to 
be applied directly at each site. 

If the coupling is nonzero, but very small, the controllability property is restored for 
n u = 2, but, according to section [V], feedback of very large magnitude is required to control 
the system due to parametric deficiency. Indeed, in order to affect the dynamics at site i 
away from i\ and 12 the control has to propagate a certain distance decaying by roughly a 
factor of e per iteration. As a result, the magnitude of the perturbation required to control 
an arbitrary site of the lattice diverges approximately as e~ Hx ^ 2 for e — > 0, resulting in the 
loss of control M. 



B. Particle in a Symmetric Potential 



The motion of a particle in a symmetric potential, such as a point charge in electric 
field, serves as another example of the relation between the structural symmetry group Q 
and the dynamical symmetry group £. This and many other interesting physical systems, 
e.g., inverted pendulum, or a satellite in orbit, are described by the second order ordinary 
differential equation 

mr = -W(r), (129) 

which can be trivially reduced to a system of first order differential equations of the form 
(0) introducing additional coordinate v = r. Suppose the potential V(r) possesses the cubic 
symmetry (described by the group O which is a subgroup of SO (3)), but is not spherically 
symmetric, for instance: 

V(r) = Vocosh(fcr) cosh(ky) cosh(A;z). (130) 
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The group Q = O defines the structural symmetry of the evolution equation ( |129| ) . Lineariz- 
ing the latter about the equilibrium position f = we obtain 



r 







r 




r 


V 











V 



(131) 



where u 2 = —Vok 2 /m, while and / are 3x3 zero and unit blocks, respectively. If Vq < 
the equilibrium is unstable, and control should be applied to keep the system close to the 
equilibrium state. 

Equation (|131| ) is spherically symmetric, with C = SO (3) and, therefore, Q C C, i.e., the 
symmetry of the linearized equation is higher than the symmetry of the original nonlinear 
evolution equation. (In fact, the full symmetry group of equation ( |131|) is C = GL(3), but 
we choose to use its subgroup CJ = SO(3), since it is physically more relevant, completely 
resolves the structure of the Jacobian matrix and, as such, correctly represents the effect of 
symmetry on the control setup.) 

Next we notice that the representation T of the group C in the six-dimensional tangent 
space {r, v} can be decomposed into a sum of two equivalent three-dimensional irreducible 
representations of SO (3) (vector representations, which coincide with the respective irre- 
ducible representation of GL(3)): 



T = 2T\ 



di = 3. 



(132) 



one needs at 



This indicates that in order to control the unstable steady state r = \ 
least three independent control parameters, n u = 3. 

Arguably the simplest way to control such a system is to re-adjust the potential (applying 
external fields, shifting support point, etc.) based on the instantaneous values of the position 
r and velocity v of the particle. This corresponds to picking the control matrix in the 
following form: 



B 




bi b 2 b 3 



(133) 



where bi, b 2 , b 3 could be chosen as any three linearly independent vectors in 



VIII. CONCLUSIONS 

Summarizing, we have determined that if the system under consideration is symmetric, 
it cannot be considered generic with respect to conventional chaos control techniques, and 
its symmetry properties should be understood prior to constructing a control scheme, even 
if the symmetry is only approximate. The failure to observe the restrictions imposed by the 
symmetry on the structure of the measured output signal will usually prevent the experi- 
mental reconstruction of the system dynamics. Similarly, an inappropriate choice of control 
parameters will result in weak controllability and, as a result, extreme sensitivity to noise, 
or even worse, complete loss of control. 

From the practical point of view, the main result of the symmetry analysis is that the 
minimal number of independent control parameters required for control, as well as the 
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minimal number of independent scalar observables required for the reconstruction of local 
dynamics, can typically be determined without any knowledge of the evolution equations 
governing the dynamics of the system. One only needs to know the symmetry properties, 
such as spatial and temporal periodicity, of the target state, and the structural symmetry 
of the dynamical equations, which in the case of extended chaotic systems is often uniquely 
defined by the geometry of the underlying physical space. One should, however, realize that 
this typical pattern does not apply to all symmetric systems without exception. The dynam- 
ical equations might, in principle, be symmetric with respect to transformations unrelated 
to "geometrical" symmetries, such as rotational, reflectional, or translational invariance. 
Additional "nonphysical" symmetries can also be introduced as a result of the linearization 
procedure. 

A number of comments have to be made regarding accidental degeneracies. We found 
that when accidental degeneracies are present, restrictions obtained using symmetry consid- 
erations alone provide only the necessary conditions for controllability. In particular, one 
obtains a lower bound on the minimal number of control parameters. Exact determination 
of that number in this case requires additional information about the structure of the Jaco- 
bian matrix, which can be gathered using experimental reconstruction. On the other hand, 
experimental reconstruction itself is only possible, if there is an adequate number of inde- 
pendent scalar observables. This number, however, is similarly undetermined. In practice, 
though, one rarely has to worry about such complications, since accidental degeneracies are 
not common and unlikely to be a problem for most actual experimental systems. Besides, 
an estimate for the minimal number of observables and control parameters can always be 
easily obtained using combinatorial arguments. Also, one should be careful in equating 
the minimal number of observables or control parameters with the highest degeneracy of 
the Jacobian matrix, especially if this degeneracy is at least partially accidental. It can be 
argued that accidental degeneracies between eigenvalues from the same irreducible invari- 
ant subspace typically will not increase the dimensionality of the respective eigenspace and, 
therefore will not lead to additional degeneracy in the local dynamics. 

We also established that it is not enough to find an adequate number of control pa- 
rameters (or observables). These control parameters (observables) have to satisfy certain 
conditions. In particular, perturbation of the control parameters should completely break 
the dynamical symmetry. The more strict independence condition is specific to each target 
trajectory and, on the one hand, requires the knowledge of the system's response to variation 
of different control parameters (which can be obtained experimentally, if necessary) , but, on 
the other hand, allows one to choose the minimal set of control parameters systematically, 
avoiding trial and error search. For example, in case of extended dynamical systems with 
local feedback the independence condition usually imposes restrictions on the mutual ar- 
rangement of "control" sites, while the number of "control" sites coincides with the number 
of control parameters. 

Finally, we discovered that the conventional approach to system identification also has 
to be modified in the presence of symmetries. In particular, in order to preserve not only 
the topology of the original attractor, but also the symmetry of the original dynamical 
equations, one has to use a number of simultaneously measured observables, which have to 
be the components of an equivariant vector function of the actual state of the system. The 
restrictions on the output can be relaxed somewhat in the case of local reconstruction in the 
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vicinity of some target trajectory. However, even then a number of independent observables 
should be used instead of just a single one, as long as the symmetry of the target state is 
nontrivial, leading to the increase in the dimension of the embedding space. Otherwise, the 
conventional approach carries over with minor modifications. 
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